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PREFACE 


Final report of Asymmetrical Booster Ascent Guidance and Control System 
Design Studies performed under Contract NAS9-13568 are contained in five 
separate volumes identified as follows: 

Volume 1 - Summary 

Volume II - SSFS Math Models - Ascent 

Volume III - Space Shuttle Vehicle SRB Actuator Failure Study 
Volume IV - Sampled Data Stability Analysis Program (SADSAP) - 
Users Guide 

Volume V - Space Shuttle Powered Explicit Guidance 
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ABSTRACT 


This report was prepared under contract NAS9- 13568 for the Guidance and 
Control Systems Branch of the Avionics Systems Engineering Division, 
Lyndon B. Johnson Space Center. The report presents an optimum 
powered explicit guidance algorithm capable of handling all Shuttle 
e-toatoospheric maneuvers. Theoretical and practical basis for the 
currently baselined Shuttle powered flight guidance equations and logic, 
as documenteo by MIT, is presented. The development is based on 
previous work by JSC, MIT, MSEC, Boeing and others. Detailed flow 
diagrams for implementing the steering computations for all Shuttle 
phases, including powered Return to Launch Site (RTLS) abort, are 
presented. Derivation of the powered RTLS algorithm is presented, as 
well as detailed flow diagrams for implementing the option. The flow 
diagrams and equations are compatible with the current MIT Shuttle 
powered flight document. 
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1.0 INTRODUCTION 

The original powered flight guidance software proposed for Shuttle consisted 
of various specialized routines to handle the various phases of flight (as 
Saturn V/Apollo). The Boeing developed Linear Tangent Guidance (LT6) 
equations (Reference 1) were baselined for Shuttle nominal ascent. Investi- 
gation revealed that Apollo type guidance was inadequate to handle some 
phases of Shuttle flight (i.e., low tnrust deorbit). It was later dis- 
covered by Long and McHenry of the Mission Planning and Analysis Division 
(MPAD) of JSC that the basic LTG algorithm was accurate and flexible 
enough to handle some of these difficult phases. Long and McHenry 
(Reference 3) then developed the concept of extending the basic ascent 
LTG equations to handle all phases of Shuttle powered flight, including 
Abort Once Around (AOA) and Return To Launch Site (RTLS) abort. Brand, 
Brown, and Higgins (Reference 2) of the Massachusetts Institute of 
Technology (MIT) later developed concepts that made it feasible to extend 
the ascent LTG to handle various phases of Shuttle flight. 

The basic Boeing guida>^ce task under this contract was to coordinate with 
MPAD, MIT, and Rockwell International (RI) in develoixnent of the Unified 
Linear Tangent Guidance (ULTG) and to implement this capability in the 
JSC Space Shuttle Functional Simulator (SSFS). Emphasis was placed on 
abort capability. Various candidate abort techniques were critiqued. It 
was determined that all of the existing candidates for RTLS abort guidance 
had certain limitations and disadvantages. As a result, Boeing developed 
a complete solution to the RTLS abort problem satisfying all constraints 
and compatible with the ULTG as a simple option (Section 5.2). In order 
to solve the RTLS problem, the original LTG algorithm was significantly 
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modified. A unified optimum guidance algorithm (References 8 and 13), 
involving higher order accuracy and more flexibility, was developed. 

This accuracy and flexibility made range throttling possible (Section 4.2) 
and made it possible for the steering to modulate to deplete excess pro- 
pellant and simultaneously satisfy velocity and position constraints 
(Section 5.2). 

The guidance work performed under this contract contributes to all phases 
of Shuttle powered flight guidance. A technique for handling the Orbiter 
Maneuvering System (QMS) ascent phase (Reference 9) was developed end is 
presented in Section 4.1. This techninue has been implemented in the 
ULT6 of the SSFS and incorporated into the currently baselined Shuttle 
powered flight guidance equations (Reference 10). Theoretical background 
and development of the unified optimum guidance algorithm are presented 
in Sections 2.0 and 3.0. Techniques for optimally handling various sets 
of end-conditions are presented in Section 4.0. (Techniques that were 
undeveloped for Saturn 'Apollo technology). Two explicit powered flight 
guidance algorithms for RTLS are presented in Section 5.0. Section 6.0 
contains a method for implementing the steering computations for all 
Shuttle phases. RTLS time-to-go and range throttling equations are 
also included. References are listed in Section 7.0, Derivation of 
generalized higher order thrust integrals is presented in Appendix A. 

A simple gravity model capable of accurately handling all Shuttle powered 
maneuvers is presented in Appendix B. However, this is not the model in 
the currently baselined equations (Reference 10). The RTLS powered flight 
option is presently being programmed in the SSFS by NASA/JSC and LEG. 
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The guidance algorithm has the following features: 


Optimality - All necessary conditions for optimality are satisfied to 
within the accuracy of the guidance algorithm. All flat earth assumptions 
are eliminated. 

Accuracy - This guidance algorithm is an approximate closed-form solution 
to the powered flight trajectory optimization problem. 

Stability - There is considerable feedbar ■» the algorithm; i.e.» every- 
thing is a function of everything else. However, in the final analysis, 
everything is a function of final state, costate, and final time. These 
final quantities converage to virtually constant values after a few cycles 
through guidance. 

Flexibility - The algorithm can handle any reasonable set of maneuver end- 
conditions in an optimum manner. The number of end-conditions can vary 
from one to seven. 

Simplicity - The basic equations are very clean and simple, e.g., cleaner 
and simpler than the Saturn Iterative Guidance Mode (I6M), the original 
LT6 approach of Reference 1 , and the equations of References 2 and 3. 

It will be seen in the following sections that the concept presented in 
this report actually represents a problem solution approach. The approach, 
although niroerical integration could be employed, eliminates expensive 
(solution time and computer storage) Newton-Rhapson type iterations. 

Desired maneuver end-conditions are expressed as functions of the terminal 
velocity and position vectors. The navigator furnishes instantaneous 
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values of state, time, thrust acceleration and gravity acceleration. 

Analytic expressions are employed to represent the magnitude of the thrust 
acceleration time history. The problem is to determine the thrust vector 
time history (and issue steering commands and throttle command if ap- 
plicable) that maximizes final mass subject to satisfying maneuver objectives 
(end-conditions). Theoretical development furnishes the optimum form of 
the thrust vector (steering parameters) and missing end-conditions necessary 
for optimality. The problem solution approach is to develop a closed form 
approximation for determining values of the steering parameters that 
satisfy all necessary and desired conditions. The solution is a predictor/ 
corrector or recursive type iteration. Several pre-thrust or pre-guidance 
cycles through the equations are assumed necessary to converge the algorithm. 

The solution approach is as follows; The solution is linearized in terms 
of steering parameters, thrust and gravity integrals and final state. The 
steering parameters are then determined from the linear equations. The 
steering parameters and thrust integrals are in turn used to predict final 
state in a manner that approximately satisfies down-range conditions. Down- 
range conditions already being satisfied, cross-range and radial conditions 
as well as appropriate magnitudes are corrected if errors exist, i.e., a 
discontinuity is introduced into the predicted state to satisfy end-conditions. 
A one pass iteration is repeated on the next periodic guidance cycle. The 
algorithm is extremely convergent, even for bad guess starting values of 
final state, costate and time. An initial estimate of velocity-to-go is 
needed to start the algorithm. 
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2.0 THEORETICAL BACKGROUND 


The guidance and trajectory optimization problem fits in the class of the problem 
of Bolza or more specifically the problem of Mayer (Reference 5, Chapter 4). This 
problem can be formulated as follov^s; consider tht class of functions 

S^(t), i = 1, n , where t is the independent variable. 

Satisfying the constraints (differential equations) 

^j(t, S. , S.) = 0, j = 1, .... t , where S h ^ , 

which involve n-^^ degrees of freedrm. These functions also satisfy the end 
condi tions 

Ck(ti, = 0, k = 1, . . . , ii 


^k^^2’ ^i2^ = 0, k = m+1, s < 2n+2. 

The problem is to find that sp'‘Cial set of these functions (S.) that minimizes 
or maximizes 


G(t, S.) 


This problem can be treated in a i imple and elegant manner by introducing a set 
of variable Lagrange multipliers, 

j = 1, .... £ (constraints) 

and a set of constant Lagrange multipliers, 

k = 1 s (end conditions), 

and employin', the principles and theorems in Reference 6, Chapter 7, and Reference 7, 
Chapter 2. 


For simplicity the notation of Reference 6 is employed here, e.g., 


n 



i=l 


j=l 
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The above problem can be reduced to minimizing or maximizing the augmented 
functional 

to 


P = [6]^ + + J Fdt, where F = Xj<}>j = 0 


In order to pxtremize P subject to constraints the first variation, 6P, must vanish. 
The following commutative lav.s from Reference 7 are used in formulating the first 
variation, 6P. 


^2 ^2 

^ L ^ L 

ti t^ 

where x is the independent variable. It is also shown in Reference 7 that 
5(F-j/F2) = (FgSF^ - as in differentiation. 

In order to obtain all of the necessary conditions for optimality shown in 
References 5 and 6, the equations are written in parametric form as in Reference 6, 
i.e. , 


t = X, S, = S.(x) from which it follows that F = F ft(x), S.(x), 

^ ^ [_ ’ dt/dx 


dS./dx 


dt , 
dx 


6F = 8F 6 s. + ^ 6 s. + 3F 6t + 3F 6t (where t * ^ dS/dx ) and 

3t 3t dx dt dt/dx 

6(^) = [(dt/dx) 6(dS/dx) - (dS/dx) 5(dt/dx)]/(dt/dx)^ 


cr 


d 

^■It 


dS 

it 


d 


6$ =-^ 6S " -^r 6t since 6 and are commutative, and 6t = -^ fit. 


With the above definitions we are now ready to formulate the first variation, fiP. 
With the usual integration by parts, removing some terms from the integral, the 
first variation is as follows: 

fiP = T + I, where 

T . [dG + (F-5 ,Fj )dt + Fj dS,]^ + U|;dC^ 

* " ' /t^ I ''5, ■ *^i * Ik '■ • 



It is necessary for optimality that the coefficients of the differentials vanish. 

The vanishing of the coefficients of the differentials in the expression T is called 
the "transversal ity condition". The coefficients of 6S. in the integral are called 
the Euler-Lagrange equations. 



Inspection shows that 

-^F - H, called the Hamiltonian ' the system. The coefficient of fit 

in the integral becomes H + Fj = 0 (called the "first integral"). Consequently, 
for problems where F is formally independent of t, the Hamiltonian is a constant 
of the system. 

H = C (constant), and the transversal ity condition reduces to 
[dG - Cdt + F^ dS^.]j + yj^dC,^ = 0. 

So, in conclusion, the "transversal ity condition", the Euler Lagrange equations, 
the "first integral", and the constraints and end conditions being satisfied form 
the necessary conditions for optimality in the problem of Mayer, stated above. 

The above necessary conditions are general and consistent with the problem statement. 



3.0 NECE.SSARY CONDITIONS APPLIED TO GUIDANCE 


For reference, the necessary conditions for optimality are summarized here in 
common notation. 



The exoatmospheric equations of motion are 


V = a^p(t) + G(R) 


(4) 


R = 


V 


m = fy 

Where is velocity, £ is radius, m is vehicle mass, and G is gravity. 


a, = X» where T is engine thrust, 
m 


P = 


sine cos ijil 

sin i|)j , a unit vector 
cose cos ^ 


(5) 


defining the thrust direction (when implementing the guidance law the small 
thrust vector control oscillations may be ignored). 


f 
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An ingenious method for handling gravity in the problem formulation has been 
developed by Jezewski (Reference 4). The approximation is made that the 
gravity vector is a linear function of the radius vector, i.e. 

G = - (p/R^) R where R is a 

constant radius magnitude chosen judiciously. It is seen that 

p/R^ = /Rp = where ui is the angular rotation rate of a circular 

orbit at radius R^. This leads to 
2 

G = - (0 R. 

It will be seen in the following formulation that this approximation eliminates 
flat earth assumptions in the guidance law. The augmented functional becomes: 

F = A •• [i - aj2 + + « *• (i " V) + e(m - f^) 

where for constant thrust m = fy = constant, and for constant acceleration 
fy = -mac/IsP where ac is the value of acceleration and I^p is engine/ 
propellant specific impulse (assumed constant in this problem), and 

{X:a:B) = (i:R.:m) h S. and (V^:£:m:e:i|») = S . , j = 1, . . . , 7 i * 1, ..., 9 

The typical powered flight trajectory and guidance problem* is to determine the 
thrust direction, P(t), time history that satisfies final conditions, 

Cj 2 (VptJ^), and simultaneously maximized final mass, m 2 - Therefore, Gem. 

Tne e^ conditions, in general, are not explicit functions of t, m, and £. 

The number of end conditions, s, can be anything from one to seven. In general, 
the initial conditions are fixed. 

Applying equation (1) 

{^) SL ~ ~ h. 

(b) a_ = h. ^ (»)^ 2l ~ ® 

(c) i = -(i-£) !!l - e 

?m am 



* The problem here is restricted to a single bum maneuver, i.e.', no long coast 
arcs between stages. 
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Referring to equation 3 


(e) H = C » constant 


(f) 


since F is formally independent of time. Therefore, the term H 
disappears from the transversal ity condition for variable initial 
and final conditions; i.e., 


Ht2 - = C - C =0 


For simplicity it is assumed here that initial conditions are fixed 
and the initial differentials vanish, therefore, the transversal ity 
condition becomes 


(g) ^2 " 


(h) I2 









since 





(end conditions are independent of t, m, e and'l'). 


Conditions c, e, and q furnish superflous information fc the problem 
as stated, since they can be satisfied without controlling anything. 


The solution to the differential equation of condition b 


(i.e. , = 0) 


is 


\ = _x(0) cos o)t + sin cot 

u 

\ - -u^(0) sin u)t + ^(0) cos wt 


In most low earth orbital maneuvers w is approximately the range angle rate. 
Since this is a homogeneous set of equations (i.e., one of the components is 
arbitrary) the above equations can be expressed as 

£ = ^ cos o)(t-K) + ^ sin u)(t-K) = J^(t) (7) 

£ = - ojX sin ii)(t-K) + i cos w(t-k) = ^(t) (8) 

where t is ijefei^enced to K » (l/ 2 )t 2 , tj » 0, ^ = ^(K) is a constant unit 
vector and \ = _x(K) is a constant rate vector. Equation 7 is the key to the 
simplicity and accuracy of this guidance algorithm, and replaces equation 4 
of Reference 1. In general the quantity ^ ~ ^ is small and |_x| » w, therefore 
_u is approximately a unit vector. "" ^ 



The thrust unit vector is expressed as 


£ * £/|£l (condition d) 


( 9 ) 


Expanding equations 7 and 9 and dropping terms of order higher than three, £ 
can be approximated as 

£ = X (1 - TX - + X (x - Tx^ - 1/6 (u|x^) (10) 

where x = t-K and T = X. • X, , T ’ = T (1 - 7/6u^X^), and = 3o)j - zJ . 

The equation of motion becomes 

( 11 ) 

The left hand side of equation 11 can be integrated twice in closed form using 
the generalized thrust integrals defined in Appendix A to yield 


-TH - 


Rjh = SpX + QpX where 
Lp = L - Tjj - 1/2 where * Oj 

Sp = S “ TQj “ 1/2 10 Q 2 » where Qj ~ Qj 
Op = Jl - TJ 2 - 1/6 CO 2 O 3 
Qp = Ql - TQ2 - 1/6 4Q3 


7/6 A 3 

- 7/6A3 


Equation 10 is probably adequate for all Shuttle maneuvers. However, if more 
accuracy is desired, a constant factor, f, can be introduced into equation 10, 
as in Reference 8. Then the average value of £ is forced to unity, i.e.: 


K 

£•£ dx = 1 
-K 



The above equation is a quadratic function of f» i.e., of the form 


af^ + bf + c = 0. 

An alternate approach is to use equations 7 and 9 and employ numerical integration 
to integrate (ay/ju|)cos wx and (ay/|uj)sin u»x. 


The right side of equation 11 can be handled as: 
^GO 

J(;i-G)dt=V^ - V - Vg . 

t. 


I (V-G)dt^ = ^ - Rg . 4, 

where ’ y and Rj are terminal values of velocity and position. Gravity integrals 
' Vg and are defined in Appendix B. The velocity and position equations 
are expressed as two linear equations in ^ and L 


r 

- n 


Lx + OA = 

()— p^ -HiN 

SA + Q.A = Erw 

nr- ^r>— 


( 12 ) 

(13) 


Solving equations 12 and 13 simultaneously and forcing ^ to unity it follows that 
i, - unit ( 4 ^ - Bg^Jp/Qp) = unit (Vjo) . (14) 

A scaler equation used to force (i.e. £j) to be consistent with TgQ 
(time-to-go) is obtained by performing the dot product of ^ with equation 13. 

Sp + QpT = where T = ^*2^ 

This constraint can be satisfied by defining as 



from v'SIch it follows that 1 ‘R'qn ® 


substituting the *'ight side of equation 15 into equation 13 (for 

i ' (Sp/Qp) C(BGN/i-^GN>-iJ * T(Rq,/X.R 5„) (16) 


♦Tl'.is same approach could be used in equation 12, possibly speeding convergence. 
First guesses for and Rgj^ could have any magnitude as long as the directions 
were somewhat reasonable. 
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Equation 16 can be bypassed if there are no position constraints. Inspection 
of equation 16 shows that is the secant of the angle between 

and X (i.e.» ^q) and that (RQ[^A*RQj»j)-x. is normal to x and has the magnitude 
of the tangent >f the same angle. So \ has a component normal to and a 
component along equal to the product of T » and the secant of the angle 
between ■^N and ho' 

We now have values of x and _x, therefore the value of Aq can be defined for the 
next guidance cycle (from equation 12) 

Time-to-go (Tqq) is determined from 



Referring to the definitions of and Equations 12 and 13, the predicted 

values of and JR^ are obtained as: 

■^RED " -g -T (18) 



■^REO "" ^p^ ^p^ ^ Bg B.J (19) 

quantities to be used in the end-condition routine in a predictor/corrector manner. 

Having determined £(t) it is possible to employ numerical integration to solve 
equation 11 yielding any desired order of accuracy for state extrapolation and 
gravity effects. This would eliminate the generalized integrals of Appendix A 
and the need of an approximate gravity model? possibly simplifying the algorithm. 
However, this is not necessarily recommended. It would possibly result in an 
increased cycle time. 


* e.g., the model of Appendix B. 


'V 


i«k > * 
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3.1 Transversality Conditions 


The following set of end conditions is defined to be used in development of 
the transversality conditions: 

Ej = - 2p/R (Energy), V = |H, R = lR|, and = V-V . 

2 2 

Eg = (RV) - (R'V) (Angular momentum magnitude) (20) 


E3 “ (M) • Uy 


(Angular momentum vector) 


Where £ is radius vector, is velocity vector, y is the gravity constant 
and Uy is the desired value of the angular momentum unit vector. These conditions 
are to be satisfied at the end of a maneuver. For convenience, final sub- 
scripts are dropped. The general transversality condition from the last 
section is: 


s 





(21) 


Where s is the number of end conditions (condition h). 
Case 1 


Consider the case where energy is the only end condition (i.e., s=l. Cj=Ej. 
Using equations 20 and 21 ^ ^ 

u = -2yV 

( 22 ) 

« 

= -2y£ 

Where 6 » - is the gravity vector. 

? 

Equation 22 implies that the final value of ^ is parallel to the velocity vector 
and the final value of is parallel to the gravity vector. 
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Case 2 


Consider the case of two end conditions, energy and angular momentum. 

s “ 2, Cj “ ^ 1 * ^2 * ^2 
Using equations 20 and 21 

"A ■ ^2 — * ^ 

‘2 - I'iIIr * Vz-^ ° 

Where = V X(RXV)and ^ = R'X(RXV). 


Performing the appropriate dot products the following scalar equations results 
-u*V - V 2 G*V= 0 
ii*£ ^2 — ^ 

-£*ij + vj ^ 

• * * 
iL*ii “ ^2 — - * ^ 

from which the two necessary conditions result 

ft 

u-V - u* 6 «= 0 (23) 

£'l^- U'Ny = 0 (24) 
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Equ.itions 23 and 24 can be written as 


k = (23) 

Sn = - (V/Ru„)u S (24) 

0 ni m a 

• • • • • 

where = Unit (u.), = [uj , = sin (u, G), 

“ liil » Up “ unit (Ny), and 
= [unit (ji)] • Up * sin (angle of attack -a). 

Inspecting equation 24 it is seen that the turning rate vector is rotated from 
the terminal gravity vector (negative radius vector) by an angle (e) that is 
proportional to the terminal anole of attack. In many maneuvers 
• 

(V/Ru_) u„ = 1, when V “ u (angular rotation rate), 

'mm R 

• 

u V 

m - -B, and u^ “ 1. In this case 
K in 


S„ “ S or 0 “ a* 
e a 

Equation 23 provides magnitude information and equation 24, direction information. 


Case 3 

Consider the case of three constraints 
s = 3, Cj = Ej, Cg = Eg, Cg * Eg 


Using equations 20 and 21 and performing the necessary algebraic manipulation, 
it is seen that equations 23 and 24 are the necessary conditions for optimality 
as in Case 2. It can also be shown that these are the transversal ity conditions 
for a velocity vector constraint (i.e., V, y* A^) with no position constraint. 

E3 can be considered an azimuth constraint rather than a position constraint. 

Case 4 

Considering the case involving position constraint(s) as well as velocity 
constraint(s), it can be shown that equation 24 is the transversal ity condi* 
tion and equation 23 is not used. 



1 
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I. 
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, r 



. ! 
j * 


In summary, three modes suffice to handle all sets of end-conditions (i,e., 
end of thrusting maneuver constraints): 

Mode 1 - Energy constraint, 

Mode 2 - Velocity vector constraints, no position constraint(s) , and 
Mode 3 - Position constraint(s). 

The transversal ity conditions are implemented in, the LTG algorithm with very 
simple modifications. Terminal values of u and ^ require several pre-guidance 
iterations as do ^GO’ -T’ and (time-to-go, terminal velocity and position 

vectors). However, after several iterations all the terminal conditions 
converge to virtually constant values (i.e., to within the accuracy of the 
guidance algorithm). The following definitions are made: 

* I-tI 

= Unit vector normal to desired orbit plane, input or calculated inflight. 
= unit (Rj) 




All of the above defined quantities are defined and available in the LTG 
algorithm. An arbitrary manner of implementing the transversal ity conditions 
is shown here. For convenience, final subscripts are dropped. From equations 
7 and 8 where t = 2K: 

jj » ^ cos wK + (^w) sin «K (7-a) 

u, * - uX, sin uK + ^ cos wK 





The various modes are handled as follows:* 


MODE 1 - Energy 
Uy « unit (yj 

it “ (ii*iiv)iiv 
it • 



*A reasonable first guess estimate for x (for all modes) results from equation 7-a 
assuming that u ■ un1t(y), I.e., f1nal~~angle of attack Is small. 
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MODE 2 - Velocity 
• • 

» unit (u) 

u = (u-^ii,,*V) u, or = unit (x) and 

. X • (6 cos u)K + Vw sin uK) , 

\ = z :::: x 

Cv cos d< - (G/oi) sin uK] 

~ “o{X Sin uK + cos <»)K 


MODES 2 and 3 - Velocity or Position 


Sg = [Unit (u)] • + u • (UyXV)/u^Rjj - 0 


iL * Ji [Ujjj U 2 - (u*U. 2 ^ii] (small angle approximation) 

and finally, i is obtained for all modes by solving equations 7-a and 7-b simul 
taneously 
• • 

X, » ^ cos uK + tojj sin o)K (25) 


The quantity X is. practically a function of K 
converge to virtually constant values after a 
algori thm. 


Tew passes tnrougn 


Equation 25 should be evaluated before and after the transversal ity computation. 
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4.0 END-CONDITIONS 


It is simple to implement this guidance algorithm to handle any reasonable set 
of end conditions in an optimum fashion. For example let the desired end 
conditions be velocity magnitude (V..), radius magnitude (Rr,)» flight path 
angle (y„)» and orbital igclination (ijj). i.e., there are four end conditions 
and the Hownrange angle and nodal angle are free. A position constraint (R^) 
is involved and it will be seen that the transversal ity condition is handlea 
as in mode 3 in paraoraph 3.1. A simple iterative manner of handling 
this problem is as follows: 

Using equations 18 and 19: 

V = VpUrU R- RpUnit 

= Unit (R), Uy = Unit (VXU^), Uj, ■■ U^KU^ 

Define an error term: 

” ®1 ^ ^D^ where a^ is a damping coefficient (a^ ± l/cos^|»)* 

empirically chosen, and n is the earth rotation unit vector. 

Then the solution proceeds as follows: 


C 

e 




Sin Yj) 



cos Y 




When 0 end conditions and the transversal ity condition are satisfied. 


Another example is in Reference 9 where the magnitude of apogee radius (Rg) is 
the only constraint. The transversal ity condition is Mode 1-energy. This 
end-condition scheme makes use of the relationship 


Rp “ [Ro(Ra-Ro)3/C(Rg/cos^Yj))-Ro]» where Rp is perigee radius.** Another 
relationship that could be used to develop end-condition schemes is 


Rp[2-(l-cose)(Rjj/R^)] - Rjj(l+cos e) * 0 , where e is true anomaly. 
* Where ^ Is terminal latitude. 

♦*Rjj and yq are computed from predicted values, Rpp^Q and 



4.1 Ascent Orbiter Maneuvering System (OMS) Burn (one constralnt- 
Apogee Control) 

A present ground rule for OMS ascent burn is to have a separate OMS guidance 
and targeting stage. During main engine burn guidance is targeted explicitly 
for Main Engine Cutoff (MECO) conditions to accommodate external tank (ET) 
impact. After MECO guidance is restarted and the OMS is targeted for 
different conditions. There is a very low acceleration during OMS burn 
( S'. 06 g's); therefore velocity direction change and altitude constraints 
could be very expensive in terms of performance and guidance and control 
(G&C) systems stability, (e.g., a velocity change of one degree at a velo- 
city of 25,000 EPS costs about 450 FPS, compared to a total nominal OMS 
propellant loading of 600 FPS. An altitude constraint could also be very 
expensive and produce erratic attitude and attitude rate cotranands). 

The guidance and targeting scheme, presented herein for Implementation in 
the SSFS, produces an optimum (minimum aV) OMS burn to achieve a desired 
target apogee altitude (or radius). Perigee altitude is not explicitly 
constrained but is implicit in terms of desired apogee altitude and OMS 
cutoff altitude and flight path angle. Predicted values of cutoff altitude 
and flight path angle are obtained by enforcing a near gravity turn 
trajectory for the OMS burn, until a velocity magnitude is attained that sat- 
isfies apogee altitude. The problem is to determine the desired velocity 
magnitude as a function of desired apogee altitude and predicted cutoff 
altitude and flight path angle. 

The three following orbital mechanics equations are used to express perigee 
radius (R ) as a function of desired apogee radius (R/^), burnout radius (Rp), 
and burnoQt flight path angle (yjj). ^ 


a “ .5 (R^ + Rp) 

(26) 

e • (R;( - Rp)/(RA * Rp) 

(27) 

cos^Yq = [a^(l-e^)/R(2a-R)3 

(28) 


Substituting equations 26 and 27 into equation 28 results in 
Then desired velocity is as follows 

*0 * ["<1; ■ alj"* O'" 
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2 

The value of cos Yh is' obtained from values of cutoff position and velocity 
vectors (Rj and V 7 . 

T 2 

.2. _ 1 _ \ . iirt'it M/ M (30) 


cosSjj = 1 - [Unit (n^) • Unit (V^)] 


Rjj is expressed as 

R[) ® lEfl •• 


(31) 


This guidance option is implemented in the following o''der: 

(1) Rq is obtained from equation 31, 

(2) cos Yq comes from equation 30, 

(3) Vq is obtained from equation 29, and 

(4) The desired velocity vector is expressed as Vr = Vp Unit (Vj) 

In Section 3.1 it is shown that a necessary condition for optimality (for the 
one constraint problem) is that the final thrust direction be parallel to the 
final velocity vector (Vj). This condition is satisfied by the following 
equatic ( (if the simplifying approximation that = 0 is made). 

X « (o)/tanu)K)[(V.j./x-Vy)-X ] 


where 
u » y/VR^ 

and V and R are initial values of velocity and radius magnitude. 

This guidance option was described in Reference 9 and implemented in the 
SSFS. The option has also been implemented in the Unified Powered Flight 
Guidance of Reference 10. 

An approximation of the above equation is 

^ s -y unit C (Rgj^/^*B(j[^ ^ ” A 3 • 

which is essentie'iy what is implemented in Reference 10. 



4.2 Range Throttling 

For Shuttle mission 3-B, it is a requirement to determine a constant thrust 
level to achieve a prescribed reference trajectory (for rendezvous pur- 
poses). Also, for Return to Launch Site (RTLS) abort, it is desirable 
to throttle to a constant thrust level to achieve a prescribed range 
from the landing site at Main Engine Cutoff (MECO). 

Equation 13 of Section 3.0 is used to develop a throttle equation. 

Spi a - R - VTgd - Rg (13) 

where R^ is vhe -Fixed, desired terminal radius vector. 

A descn’ption of the method is simplified by making the assumption that 
T = = 0. Performing the vector scalar product of with equation 13 

results in 

S “ - -^N • 


In order to satisfy equation 32, it is necessary to know the proper value 
of thrust, mass flow rate, and T^q (time-to-go). Thrust (F) and mass 
flow rate (m) are constant and directly proportional to each other, i.e., 

F/m = V (constant exhaust velocity). 

TgQ is inversely proportional to thrust, or visa versa. 

The authors of Reference 10 developed the concept of analytically 
determining the time-to-go (T^q) to satisfy a range equation (similar to 
equation 32). They then used^^this new value of TgQ (along with the old 
value, TgA) to define a proportionality constant (K7 relating desired 
values ot'^F and m to prescribed reference values and mj^), i.e. 

F = KFj^ and m = Km^. (33) 

The value of K is initialized es K' = 1.00 and on subsequent iterations 
is altered as 

K = K'T^/Tg and k' is reset as K‘ » K, (34) 


where T^ is the burn time until an acceleration limit is attained, or 


I ^ 1 Q 

Tg * Tqq if an acceleration limit will not be attained. 


Tg is the value of 
Tg on the previous iteration and Tr = To + ^^60 

ired increment in to satisfy the 


is the desired increment in Tg« to satisfy 
range equation. The value oOTgO is obtained from a Newton-Rhapson 
iteration employed in satisfying the range equation. Each guidance 
pass represents one iteration, however, the solution is virtually converged 
after a few passes. 
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The approach described in the above paragraph is taken here; however, the 
solution for aTqo to satisfy the range constraint is approached in a 
different manner. A value of is found to satisfy 

equation 32, i.e.. “ ® “ 

^p ~ — ^ " (a* » (35) 

where e is an error to be driven to zero. 

For small values of aTqq* a very good approximation for Sp is 

Sp ’ s; ^ .5L (Tjo-T^o) where is the uncorrected value. 

Another good assumption for variable thrust is 
3L 

= 0, since is essentially an average acceleration multiplied 

^'60 P 

by TgQ, and if the acceleration is changed by a given factor, T^q is 
^ changed inversely by the same factor. 

Therefore a good approximation is 

as. 


3T 


P 


GO 


Another approximation used here is 


3R 


(35) 


“ )L» where R and V are gravity integrals. 

-9 

Equation 35 is expressed as 

E * Sp - L*(Rj-R"VTQQ-Rg) . 

Differentiation of equation 35 yields 

If — " • 

^'go P - 

Adding to equation 36 and substracting the same quantity yields 
3e 


(36) 


3T, 


GO 


• -SL. - 1-(Vt-V-V.) t l-lf 
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Using the identity * 


Ij ■ V. “Ig ^ ■ Lp, it follows that 

|s- = X-V. - .5L . (36) 

^'go P 

And finally, 

algQ" * 


Where £ is defined by equation 35 and -f — » by equation 36. 

'"go 

Equation 34 is used to define K {where " TA + ^Tg^) and the desired 
value of F and m are obtained from 
equation 33. 

The only modifications necessary to implement this method into the range 
throttling routine of Reference 10 is to make the following definitions. 

goz p 


Ar 


Atgo 


= 1-Vt - 


•®''go2 ' 


- R - n^-Rg) - Sp . 



5.0 RETURN TO LAUNCH SITE ABORT (PROPELLANT DEPLETION CONSTRAINT) 

The objective of Return to Launch Site (RTLS) abort guidance is to provide 
closed loop pitch and yaw commands to achieve explicit burnout conditions 
at External Tank (ET) fuel depletion. The explicit conditions at Main 
Engine Cutoff (MECO) are relative velocity (Vp), relative flight path 
angle (yr), altitude (h), relative flight azimuth (A2r), range from 
landing site (Rls) burnout v/eight (Wgo)* A 36 acceleration limit is 
also imposed. The RTLS targeting algorithm is not presented in this report 
since it is presented as one of the targeting options in Figure 5-11 
of Reference 10. Equations to constrain vehicle attitude at MECO are 
easily developed, should this become a requirement. 

This Section presents two powered flight guidance algorithms capable of 
handling the burn to MECO portion of Shuttle RTLS abort. Throttling 
capability is provided to constrain range from the landing site at MECO. 
However, throttling is not necessary should it be desirable to relax 
the range constraint. The work presented In Section 5.1 is based on the 
concept developed by Russ Sievers (Reference 11), called Launch Abort 
Guidance Simulation (LAGS). LAGS is presently being used in the Space 
Vehicle Dynamic Simulator (SVDS) for RTLS performance and targeting 
studies, and is currently being programmed in the SSFS, 

The LAGS guidance concept is to have a constant attitude maneuver before 
the nominal guidance maneuver, and the value of the constant attitude is 
computed explicitly, simultaneously with the parameters of the nominal 
guidance phase (i.e., the constant attitude phase is not open-loop). 

The duration time of the constant attitude phase is proportional to the 
excess velocity or excess propellant. Time-to-go is an explicit function 
of desired burnout weight. 

The guidance algorithm presented in Section 5.2 requires no constant 
attitude maneuver to expend excess propellant. The algorithm represents 
an explicit solution to the seven constraint guidance problem (with 
constant thrust level before acceleration limit)— burnout weight is ex- 
plicitly constrained as well as terminal position and velocity vectors. 
Burnout weight is controlled by bum time, the steering modulates to 
expend excess propellant (while satisfying radial and crossrange position 
and velocity constraints), and a constant thrust level (before acceleration 
limit) is detemiined that satisfies the downrange constraint. 

The "single maneuver" algorithm of Section 5.2 has advantages over the 
"two maneuver" algorithm of Section 5.1. It is much simpler to implement 
as a guidance option in the nominal guidance (with the addition of a 
negligible amount of equations and logic). It also eliminates most of 
the disadvantages of other candidate powered flight RTLS guidance schemes 
(a comparison of candidate schemes is out of the scope of this report). 

Atmospheric terms are not explicitly included; however, an empirical 
velocity bias term (that ranges from full value to zero at MECO) is intro- 
duced to insure that enough propellant is available to achieve desired 
velocity and burnout weight simultaneously. This bias term also compensates 
for other approximations in the guidance algorithm. Time-to-go (Tqq) 
computed explicitly as a function of nominal desired burnout weight. 

If more propellant is available than required to achieve the velt-ity, 
the steering modulates to deplete the excess propellant. 
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5.1 Derivation of Constant Attitude Equations 


The RTLS guidance equations will be developed here in the order of 
computation. The following definitions are made: 

V - Engine exhaust velocity 

V- A 

a Measured acceleration 

- Acceleration limit (e.g.» 3 g's) 

W - Current vehicle weight (onboard estimation) 

Wgo - Nominal MECO burnout weight (input) 

Time-to-go is computed as follows: 

Tp = T(l-a/aj^) - (Time of constant thrust) 

( \s “ (Time of constant acceleration) 

Tqo '•‘p + (Time-to-go) 

Additional parameters are computed as follows: 

Velocity available: 

V, = [T/(,-Tp)] t 

where L„ is the velocity increment during the reorientation maneuver 
(computed later in the order of computation). L* is initialized as zero 

u 

on the first pass. Velocity required is updated as: 

*ft£Q = %Q - Sensed 'iisensed velocity Increment 

since the last guidance pass. The excess velocity Is now computed. 

l-c * Vj^ - iVj^^ql and the time of the constant attitude maneuver is 
computed as Tj « t [1 - exp (-L^./Vgj^)]. 

( 


26 



The equations of motion assume the following form: 


-1 “ ®1-1 ^1 

= 32 ! [1 - l/2u^(t-K)^]A + [(t-K) - l/6u)^(t-K)^]x ( + G 

0 1 t T 2 where T 2 = Tg^ - . 

Where 

a-, = (F/m)^ 
ag = (F/m)2 

F - Thrust (function of time) 

m - Vehicle mass (function of time) 

iC^ - Unit vector defining constant attitude phase 

G - Gravity (function of radius vector) 

2^ - Unit vector defining mid-point of nominal guidance phase 

- Pitch rate vector during nominal guidance phase 

For simplicity the following constraints are introduced into the above 
equations 

X • X = 0 
2 

u = 1 * 2 l 

The expression for was derived in Section 3.0 (assuming that 0). 
The following integrals are defined: 

/i 

j a^dt ® (excess velocity) 

0 

= (T^ “ T )L^ + VgJ, 

0 0 


Where a^ = 






' < 


t 


f 

\ •- 


i 

I 



.. » 

^ i 


The following integrals are defined in Appendix A. 
T 

L- /a^dt 

0 

T t 

S-J/ajdt dt 
0 0 

T 

J„, = f a2(t-Kfdt 

0 

T t 

%i‘I dt, m=l, .... 3 

0 0 

Where tg = t-T^ and 82 = Vg^^T^-t) or a2 = aj_ 
V2(VT,). 


Gravity integrals are defined in References 2 and 3 and Appendix B. 




% 


.fut 

0 ^ 

=/ Jidt dt. 
0 0 


Using the above integrals the equations of motion are integrated twice to yield 

LcC, * Lpi ■ - i - y, - dpi a 4, 

Sc£i + SpX + Qpi + LpC,Tj • B:f - R - VTgp - Rg a Rg^ 


Where 

Lp = L - 1/2 A 2 
Sp T S - 1/2 Jq2 

Jp * 0^ “ 1/6 u 

Qp - q, - 1/6 J Qj 


V 
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Vj is terminal velocity vector (at MECO) 
Ry is terminal radius vector 


R is current radius vector 
^ is current velocity vector 
The integral can be redefined as 

^c ' ^^1 ” ^^**c * '^ex^l ^c^2 * ^^6o ‘ ^^^c '^ex^l 


The above equations now assume the form 




S^Ci + SpX + QpX-f^ 




(38) 

(39) 


The unknown scalar factor f is introduced because it Is assumed that the 
direction of R^j^ is known but the magnitude is an unknown function of the 
thrust integrals. The desired terminal radius magnitude (R^) is an input, 
however, the down range component of the terminal radius vector (Ry) is 
unknown. Even if throttling is used to constrain down range position, it is 
assumed to be unknown in solving equations 38 and 39. At another point in 
the equations, the thrust level is corrected to correct the difference 
between desired radius vector and predicted radius vector (based on 
equation 39). 


Equations 38 and 39, in component form, represent six linear equations In 
ten unknowns (C<j, y it and f). However, four scalar equations are 
Introduced, making the solution possible in a recursive manner. Equation 38 
Is solved Independent of equation 39 for C<| and Then values of C;| and ^ 
are substituted Into equation 39 from which ^1s solved. Performing the 
appropriate vector dot products the following scalar equations result. 


Li * 




pi- - i^l 


Lei * ^1 * Lp 


Ci^L^p.Ve^.VeNHV 


* ^1 ’ 


2 

GN 


" i * Vqn ' ''gn^x 


Sgi • C, + Sp • fv • 

since C| • Cj ■ 1 , X. • X, « 1 , and ^ i ■ 0. 
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Defining 1 * £j '•t follows that: 


'a ■ ' - tC-c * ■ ''VVp 

(40) 

Cp ■ C-c * V.'/'gn 

(41) 

'x ■ '^c^a ♦ S>''''SN 

(42) 

t ■ (S,C^ + Sp)/(X • Rg^) 

(43) 


A solution to equation 38 Is constructed by making the following assumptions: 
C|» X, and ^ all He In the plane defined by 

Uy « Unit (Vgnj) and 
(L- Unit [(Bsn/!!* • Bgn> ■ 


Another essumptlon Is thet • N > 0 and x • N < 0. It 1s seen that C 

I ““ — . — a 

is the cosine of the angle between ^and Is the cosine of the angle 

between and ^N’ and Is the cosin'? of the angle between >^and hr 

The value of Is determined from equation 40 and substituted Into 
equations 41 and 42 for C. and C. . Then the solution for C, and x follows 

V A 

as: 

c, - M 
X ■ - ( \X^)N 

The steering command during the constant attitude phase Is 

(45) 

Now that X Is known the value of f Is determined from equation 43. 

In solving equation 39 for x. It Is desirable to force the radial magnitude 
constraint, {R j| ■ R^. This Is done In the following manner: 




)N 


iAA \ 


O 
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- R + + R^ 

-C — — Go 

- R„ umt ((Rg„ + R^) 


Sgn ■ % - ^ 


f (ScC„ + Sp)/(x-R|,„) 

The solution for \ (from equation 39) follows as 

i '(fBGN - SA - Spi)/«p 

The following parameters are now computed for the next guidance pass: 


O 


C' 


(Initial gqidance command of nominal guidance) 

^ coswK - (^w)slnwK 
c_ » L, 

os -0 —I 

Kunit vector halfway between and C^j.) 

£j • (£„ + c,)/>/2(i + 

(Duration time of reorientation maneuver) 

Tg • [ 1/2 - s1n"’ 

Where td|^ Is pitch rate limit (e.g. 10^/sec) 

(Velocity increment during reorientation phase) 

Le • -T,)/(t -T, - Tg)] 

(Velocity loss factor due to turning -approximate) 
d - 1 - l/6«,\ (y2)2 

The expression for d is exact assuming constant average acceleration (a * L.^/T^) 

Tg 

dig • f a[l - l/2«,\(t - Te/2)^]dt, 

The velocity required Is now computed for the next guidance pass. 

VGo-V^-V-Vg.JpiML-Lp)! 

^EQ * ^0 ■ '*■ i^BIAS 

2 

where ^ ~ ® •" empirical bias and 

T^ Is a reference time. 
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5.2 Derivation of Steering Modulation Equations (for fixed burnout weight) 

In the following analysis, Tqq and the thrust intearals L, S, 

(i » 1,..., 3) are fixed (for a fixed thrust level), since the burnout 
weight is fixed. Therefore, the steering must modulate to burn up excess 
propellant, as well as satisfy position and velocity constraints. It is 
assumed that excess propellant is available, whether the amount is con- 
siderably large or approaches zero (the amount converges to virtually 
zero at thrust termination). 

The following equations are developed in a second order manner (rather 
than third order) in order to be compatible with the equations of Reference 10, 
which are the currently baselined Shuttle powered flight guidance equations.* 
The author feels that second order accuracy is adequate to perform the 
RTLS powered flight maneuver. 


These equations also supercede the equations presented in Section 3.0 of 
Reference 12. In Reference 12 it was assumed that T = ^*3^ “ 0. However, 
since the writing of Reference 12, the author feels that the steering rate 
must modulate in direction as well as magnitude (i.e., T;< 0). The value of 
T is determined to satisfy velocity and position simultaneously. 


Dropping third order terms, the steering equations from Section 3.0 
(equations 12 and 13) are written as 


+ Jni ' Vtm 

P^ P^ 


(46) 


S^ + 
P- 


Q ^ 


hn 


(47) 


where 


hn’ ^lAS ^^lAS defined on page 31 ) 


L • I - TJj - (48) 

Sp - S - TQj - .Si\ (49) 

Op ■ J, - TJj (50) 

Qp . Qj - TQj (51) 

•2 • • 

where x*” ■ . 


The above equations are simplified by the original LT6 constraint 
Jj • 0 - U - 0. or K • 0/L . 

The above relationship is defined on page A-4 of Appendix A. Equations 48 
and 50 become 

* The f(«rerunner (presented in Reference 14) of the current steering law 
was developed as a second order algorithm. 
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X 


L, « L - .5x‘J9 
P 2 


Jp = -TJ2 


With the proper substitution equations, 46 and 47 become 
(L - .5x^J2) (-TO 2 ) i= VgN 

(S - TQj - .5x\) X + (Qj - TQg) X » fRg^. (47) 

where f is a scalar quantity to enforce downrange position. 

Performing the vector scalar product of x^ with equations 46 and 47 it 
follows that 

x^ + 2T^ = 2 (L-x.y^.J/J, = J (52) 


^ - '“p 


S - .5(x^+2T^)Q2iS - .5 u)|^Q 2 = fX-Rg,^ 


The value ofup is a well defined quantity and is the first RTLS steering 
computation, i.e., 

«P = 2(L-x*Vg^)/J2 . (54) 

The present estimates of and ^ are available and Rgj^ is updated as 
f ® (S - •5upQ2)/(l‘Rgj^) (from equation 53) 

Bgn ■ fBsN (55) 


Equation 55 provides a stable and accurate prediction of R^», for a 
fixed thrust level. 

Two relationships come from equation 52: assuming that the value of T is 
known, 

= u.^ - 2T^, and (56) 

P ,2 

assuming that the value of a is known, 

T » ± C.5(u^ - (57) 
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If both T and are known, the values of L , S , 0 , and Q are determined 
from equations 48 through 51. P P P P 

Now, equations 46 and 47 (incorporating equation 55) are solved simul- 


taneously for x^ and i.e., 

X = Unit - R^N'^p/Qp). and (58) 

i HRqn - SpX^/^P 

So, the value of x'^from the position equation (equation 47) is 

X^ » x*i . (60) 

and the desired unit thrust vector is 

= Unit C(l-Tx - .5xV)^+ (x-Tx^)xJ, (61) 


where x » At-K, and At is the time increment since the present values of 
\ and ^ were computed. 

ft 

It is obvious that to have a solution the values of x^ from ecuations 56 and 
60 must be. equal , i.e. , 



This relationship can be forced in an iterative manner, e.g.; 

X^ = .5(u^ - 2T^ + Xp) (62) 

ft 

where Xp » x^«x^ (equation 60). 

Given this value of x^, the value of T can be determined as 

T = -[.5(up - X^)3 (equation 57) . (63) 

It is obvious from equation 57 that two solutions exist to this problem: 
a positive value of T, and a negative value. The minus sign in 
equation 63 is arbitrarily chosen. 

Values of T and Xp must be assigned in order to start the algorithm. 
Arbitrary values are; 

T - 0 and x^ » . (64) 

P P 
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Using the above equations, the solution can proceed in the following order: 
• 2 

1. (equation 54). 

2 2 

2. If the first RTLS guidance pass, T » 0, and x‘ = 

3. A (equation 62) 

4. T (equation 63) 

5. S , J , and Q (equations 49, 50, and 51) 

r H r 

6. (equation 55) 

7. X (equation 58) 

8. ^ (equation 59) 

9. r (equation 60) 

10. (equation 61) 

11. Xp = x^ (for next guidance pass) 

If the steering vector routine is implemented in the proper manner, 
steps 1 through 4 are the only additional steps required for powered 
flight RTLS capability (as in Section 6.0). 



6.0 GUIDANCE IMPLEMENTATION 

In this section, emphasis is placed on implementation of the Return to Launch 
Site (RTLS) "single maneuver" abort option. Repeated reference will be made 
to Reference 10, the currently baselined guidance equations, since the RTLS 
equations presented here are tailored to be compatible with the Reference 10 
equations. The intent is to implement this RTLS option in the SSFS deck in 
which Rockwell incorporated the Reference 10 equations. 


The nomenclature used in Reference 10 
is defined here; 

\q " ^GO 
r = R 


V = V 



(different than that in this report) 

r = R 
gray g 

''gray = \ 

''gn * ''gN 

P=Q2 

H = Jg 

Q = Qi 

-TH = Jp 

L - -5X^H = Lp 


The various maneuver modes are defined in Reference 10: 

^roode “ ^ Standard ascent 

^mode * ^ R®^^®**®hce trajectory - range throttling 

^mode * ^ 6»‘«v1ty turn/apsis control (OMS ascent) 

^mode “ ® R*^urn to Launch Site (RTLS) - range throttling 





The changes to the Reference 10 equations necessary to implement the RTLS 
option will be listed block by block referring to the flow diagrams of 
Section 5 (Reference 10). 

Block 1 Initialization 

(1) After r^^gy is defined, 

V = 2r 
-qrav -grav 

T = 0 

X = 0 

(2) In = 5 initialization (RTLS) 

after is defined, 

i^ = unit (r^) 

Do Block 3.521 v^equation 
*<i = id ^ Vid 




Block 3.1 Time-To-Go 




On the line between (Do until i > n-1) and (If n > 1), 
If = 5, do the following computations: 

■'l ° ’k ' ''ex.k''^ 


If T,< 0, T, = 0 



THEN 

V ' ''«,k'" t'k/l’k - 

If T,> Tj 






‘'l ° ''ex.k’" l^'k/l'k • 

! 

ELSE 
L— - 1 

'■''2 “ ''ex.k’" <"’k9o"T.k/*'B0\> 



''go = '"'l " '"'2 


''bias ' ^ ■ ^ 0 ^^^ 


WgQ is input desired burnout v^eight. 

is reference gravitational acceleration. 

B is an input empirical bias (possibly, B=0) 
t^ is time RTLS is initiated. 


After K = J/L 
P = P - 2QK + SK^ 
H = H - JK 
Q = Q - SK 


Block 3.2 Integrals 





Block 3.3 Turning Rate 


. unit (^q) 

2 *2 

1,2,3, 

4, X ■ u/vr^ 

■ — ‘ — ^go " ^rav ^ ^go^^ go^ 

^ " ^rav ^^go^^ go^ 

igo ' '•i 

Xga ’ V ^ ’bias 

■ 6. _i> Jp = 2(L - X. • Vg„)/H 

Ipass ' ‘ 

*pass " ^ 

X = [.5 (»^p - 2T^ + x^)]’^^ 

If x*' > »^p. X . 

T = -C .5 (»^ - x^)]’''^ 


2 

Note: Up is a quantity, 

not the square of 
a quantity. 


; - TQ - .Sx'^P 
) - TP 

)da ' “ • Igo ' V ’ ^ "p 
,de ' • igo = 

. 0 . Qp = 10® 

,de ° ® 'V * %o "/%> 

Cgo - Spi)/Qp 


unit [(1 + TK - .5 xV) X - (K + TK^)i;i 





1 


Block 3.4 Predictor 


^thrust ' - 


.5 X^H) X - THX 


I I' i 


f 


thrust “ ■^^o 


Eliminate the computation. 


Block 3.53 Range Throttlinc 


= L - .5 
goz p 


If = 2 or 5 (i.e., only one computation), 

'^'■itgo ' i • id - •= “goz 

‘fz = i • (id ■ I - i‘go ■ ^>-8v> • 


I V 

f .-'f «•. 
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APPENDIX A 


GENERALIZED THRUST INTEGRALS 


Derivation and implementation of the LTG guidance concept requires the definiti 
and evaluation of the following integrals: 



m**0f If •••§ 3 


1 “ f « • • f N^ 

Where m denotes the order of the integral, is 'the number of the first 
thrust phase and Mj is the number of the last thrust phase. 

a^ = (F/M)i = Thrust acceleration of phase i 

Tq = Initial time of phase i 

Tp * Final time of phase i 

To = Bum time of phase i (To = Tp -T^ ) 

Bi F^ 0^ 

t and S are variables of integration 
K is the time about which the LTG expansion is made. 

(K is approximately the midpoint of the maneuver.) 

Evaluation of K will be shown later in this appendix. 



For the Space Shuttle vehicle it is assumed that either (thrust acceleration) 
is constant or (thrust) is constant* therefore, these integrals can be evaluated 
in closed form. After the above integrals are evaluated the total integrals of 
phase i are evaluated as: 

“ Jp ni* 0, 1* •••* 3 

mi 

i-1 

= 2 ■'no 

j--Ks 

Then the total integrals for the maneuver are: 

Ns 

•^m * 2! **** ^ 

i=Ks 

Ns 

^m “ ^ %)i ' 

This results in a total of eight integrals. 

For constant thrust, the acceleration of phase i is of the form:* 


a 


1 


Ve. 


T.-t 


Where t = S-T. 




F = Constant thrust 


m^ = Initial mass 
0 


m a Constant mass flow rate 


V- * Engine exhaust velocity 
1 

^ i ” 

*For convenience, the subscript "i" is dropped on F, m^ and A. 

A-2 



o 


The first constant thrust integral is; 


P . 
mi 




i * 0, 1 , , 3 


dt 


The following identity makes all the integrals recursive in terms of 0^ and Q„ 

^oi ^oi 






j J dt ’ I 1 

• T.-t 0 


m-l\ 


dt- 


Ti-t 


Ve^(t+T^-K) 
1 


m 


n I 


m 


Therefore » 


^mi ^ °i P(m-l)i 


- Ve, 


- (Y)”' 


m 


(1) 


C 




m 


(Tb +T - (T 

«i 0^ 0^ 




>(T- -K)\ 
°i ^i 


( 2 ) 


o 


The integrals are recursive in terms of: 

T„ 


'B 


i 


dt a Vp.ln/ 'Ml and 


0 



r t 

r Ve.dS 

0 

i vT" 

L 0 ’ J 




since (t+T^^^-K)"’’® « 1. 


i 


dt » (Tg -xjLp +V Tg . 
“i ^ M ®i “i 
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o 






For constant acceleration a. = a, = constant and 




\ (t+Tjj -K)"’dt 


"ml 


Jp = a, 
^mi ‘•i 


(Tg +T - (T 

b. 0^ 0^ 


nrt-1 


Qp = 

mi m 




IhT 


tn*“ Oj If ••ftf 3^ 

'•''i : \ ■ W 


Sd r Qp , “ ” •SLjj Tj 


i - "01 


i “i 


Pi'Bi 


K is evaluated from the LT6 constraint 
J-LK 


\ 


From above equations 
Ns 

i»Ks i-Ks 0 

Ik, Ni, 

■ 1\ ■ : '-'''• 
1’Ks I^Kj 


from which K = J/L. (Jp^ is defined on next page). 


(3) 

(4) 
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For constant thrust 

(t+T^ ) = (t^+T^ )*-d 


/ 


0 


't- V,t 


and for constant acceleration 


B.. 


Pi J ‘-i ®i h ®i h °i “i 


0 

= Lp (.5Tg +T ) 
Pi “i °i 

Where Ln = a, To 

Pi H “i. 


K is evaluated as above before the remaining integrals are evaluated, since 
K is needed in the integrals w... . m > 0. 


From inspection of equations 1, 2, 3 and 4, it is seen that the constant thrust 
integrals and the constant acceleration integrals can be put in the same form 
by making the following definitions. 


If Constant Thrust 

If Constant Acceleration 

E = 0 

E • 1 

F ® T «K 

i 0^ 

d 

O 

A » -V^ 
®i 

A » a, 
*-i 


’’ 1 = 


P'j. 




K) 

nH-E 


m+E 


In the following equations, the m+1 subscript is used to avoid the zero subscript. 


for computational purposes (e.g.» Jo s On 

Poi Pli 


3i 


= Op ) . 
4i 


9 • • • » 




( 5 ) 


p* -p* 

Up » FJp + A(-l^) 


P 3 “ P‘2 
Pi “ P’l 



(nH-l)i 



-K) 




in* If •••f 3 


The follov/ing integrals are also used 



( 6 ) 



APPENDIX B 
GRAVITY INTEGRALS 


A sipiple, untested gravity model is presented here. The model is probably adequate 
for all ShuUle maneuvers, and could have advantages over existing models. 

Assume that a quantity (e.g., radius vector) can adequately be represented as a 
third order polynomial in time. 


R(t) = R + Vt + 1/2 + 1/6 Vt^* 

(1) 

Differentiation yields 


V(t) = V + Vt + 1/2 Vt^ 

(2) 


Initial conditions R and V are known as well as final conditions T^q, Rj and V^. 

Substituting final condition into equations 1 and 2, coefficients )}_ and are 
readily obtained. 

Using these coefficients and defining 


h 




R(t)dtdt • Ry = _i_ 
^GO 


/ 


GO 


R(t)dt 


it follows that 

Ry » R + 1/2 V Tqq + 1/2 - R - yjgQ) - 1/12 (V.^ - V) Tqq 


and 

R^ = R + 1/3 V TgQ + 3/10 (Rf - i - ^60^ * ^^T " “^GO. 

Rearranging these equations results In 

^ - 1/2 (Rj + R) - 1/12 (V.^ - V) TgQ (3) 

R^ «= 1/10 (3R^ + 7R) - 1/15 (V^ - 3/2 V) Tgg (4) 

These are mean values of R to be used in gravity computations since gravity can be 
approximated as: 


* For convenience, parentheses are omitted on fractions in the follov/ing 
equations in this section, e.g., 3/2V = (3/2)V. 


are cun«;tants and 




6(t) 


£(t), where p and R 


^GO 


^GO ^ 


/ uu o « VIU ^ 

G(t)dt and 1/2 ^TgJ = J ^ 6(t)dtdt . 

From the above equations it is seen that 


^ = _jj_ Ry and ^ = -R- Br SOO'l approximations and gravity effects 

rJ 'r 


are expressed as 

V, = VgO %io 




i ^ 

! 

i'i 

i’ ' 


1 


i 

i ’ 




Hiwever, simulations were conducted by the Guidance and Dynamics Braitch of MPAD* 
using the above model. It was determined that deriving Vg and ^ in cartesian 
coordinates as above is not adequate for long, low thrusPburn ^cs such as 
Shuttle de-orbit. 

The above model is quite accurate for short, high thrust burn arcs; however, 
the "generalized" g.avity model presented here makes use of polar coordinates 
in order to more adequately handle long burn arcs. It could be thought of as 
a third order "conic + thrust" solution. Using the method described above, the 
gravity model is as follows. 


Transforming initial and final state to polar coordinates 


£ = 


% 


sin“^ K-Uy) 
sin“^ (4*U2 !j 


V-R 

1-h 


i = 
R 


^ unit (R) 


, where ^ * unit (i^^Uy) 


D 


• 

’vt-Rt" 

1_ 5 

Q 

= 

n 

. V 

0 

tt 


J^T 




Ju 


Expressing the mean values in polar coordinates (equations 3 and 4): 
fiv ' - 1/12 (ix-i)TQo. and 

» I/IO (3£^+ ?£) - 1/15 - 3/2£)Tj. 
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‘go 
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Transforming polar mean values to cartesiari coordinates. 






Ry = [cos Py COS Py + sio Py ) Uy sin Py^] Py^ , and 

2 3 3 


'2 


1 


5r ' [“= PRj (u^ cos s)n p„^) + sin P„^] 

Gravity losses are now expressed as follows: 


iv = ^ (^v> 

’ ^ 

Sr ” f (R|)) 

■ 


'GJ 


wiiere f (JR) is an accurate gravity model or e simple spherical gravity 
model (i.e., 6 = - R). It is possible that an accurate gravity model is 

R^ 

more adequate for low thrust (QMS) burns (for accurate state extrapolation). 
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